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Hadron properties in dense medium are treated in a unified way in a skyrmion model constructed 
with an effective Lagrangian, in which the p and uj vector mesons are introduced as hidden gauge 
bosons, valid up to 0{p ) terms in chiral expansion including the homogeneous Wess-Zumino terms. 
All the low energy constants of the Lagrangian - apart from the pion decay constant and the vector 
meson mass - are fixed by the master formula derived from the relation between the 5-dimensional 
holographic QCD and the 4-dimensional hidden local symmetry Lagrangian. This Lagrangian allows 
one to pin down the density rij /j at which the skyrmions in medium fractionize into half-skyrmions, 
bringing in a drastic change in the equation of state of dense baryonic matter. We find that the C/(l) 
field that figures in the Chern-Simons term in the 5-D holographic QCD action or equivalently the 
to field in the homogeneous Wess-Zumino term in the dimensionally reduced hidden local symmetry 
action plays a crucial role in the half-skyrmion phase. The importance of the to degree of freedom 
may be connected to what happens in the instanton structure of elementary baryon noticed in 
holographic QCD. The most striking and intriguing in what is found in the model is that the pion 
decay constant that smoothly drops with increasing density in the skyrmion phase stops decreasing 
at ri]^ /2 and remains nearly constant in the half-skyrmion phase. In accordance with the large Nc 
consideration, the baryon mass also stays non-scaling in the half-skyrmion phase. This feature which 
is reminiscent of the parity-doublet baryon model with a chirally invariant mass mg is supported 
by the nuclear effective field theory with the parameters of the Lagrangian scaling modified at 
the skyrmion-half-skyrmion phase transition. It also matches with one-loop renormalization group 
analysis based on hidden local symmetry. A link between a non-vanishing tuq and the origin of 
nucleon mass distinctive from dynamically generated mass is suggested. We briefiy discuss the 
possible consequences of the topology change found in this paper on the forthcoming experiments 
at the rare isotope beam machines under construction. 

PACS numbers: 12.39.Dc , 12.39.Fc, 21.65.-f, 21.65.Jk, 



I. INTRODUCTION 

Understanding the QCD phase structure at high den- 
sity is a wide open problem in nuclear and particle physics 
as it involves the highly nonperturbative region of strong 
interactions so far unaccessible by lattice gauge calcula- 
tions. Any theoretical investigation of the QCD phase 
structure, whether it is at the level of models or effective 
field theories or more fundamental nature, would there- 
fore be highly desirable. In this paper, we investigate 
what happens to hadrons at high density in the theoret- 
ical framework developed in Refs. [llQ. To access dense 
baryonic matter, we exploit the approach proposed in 
Ref. [3| where both the elementary baryons and multi- 
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baryon systems are described in a unified way with a sin- 
gle effective mesonic Lagrangian with appropriate sym- 
metries of QCD, the most prominent of which is chiral 
symmetry, from which baryons are generated as solitons 
(called skyrmions) and many-body interactions are taken 
into account with the given chiral Lagrangian. 

In the past, for example in Ref. [J], dense matter was 
studied with the Skyrmc model put on crystal lattice 
with pion fields only [5| or minimally implemented with 
the lowest-lying vector mesons p and u that will be re- 
ferred to in what follows as HLSminl^r, p, tj)lj. There are 
however no reasons to think that such models would suc- 
ceed to describe adequately, not to mention quantita- 
tively, dense baryonic matter. The Skyrnie model has 
the leading current algebra term which can account for 
very low-energy dynamics but the Skyrme quartic term 
that is crucial in stabilizing the soliton is neither unique 
nor sufficient. Given other degrees of freedom that can 



This terminology will be defined precisely below. 



enter, there can be a variety of different ways of stabiliz- 
ing the sohton. Even at the leading order in Nc, there is 
no reason why there should not be other quartic as well as 
higher derivative terms. But accounting for them - which 
in principle should be feasible - will bring in a plethora of 
uncontrolled parameters, rendering the approach totally 
devoid of predictive power. As for the minimal vector 
meson-implemented Skyrmc model HLSmin(7r,p, w), the 
approximation used to simplify a part of the hidden local 
symmetry Lagrangian, namely the homogeneous Wess- 
Zumino term (hWZ term for short), is unjustified when 
the vector meson mass could decrease in density, which 
is the basic premise of hidden local symmetry (HLS for 
short) theory ^^. 

In this paper, we approach the problem with the 5-D 
Yang-Mills action in curved space plus the Chern-Simons 
term that has been arrived at both dimensionally de- 
constructed, bottom-up, starting from the current alge- 
bra term Q and, top-down, from gravity-gauge duality 
in string theory [l^, Ull D Kaluza-Klein- reduced to four 
dimensions, this theory consists of an infinite tower of 
hidden local symmetric gauge fields. Our assumption is 
that this is a theory of QCD type valid in the large Nc 
- and in the case of the top-down model, large 't Hooft 
A = g^Nc - limit. One may ultimately use the 5-D action 
for dense baryonic matter as was done in Refs. |12| - [l4j 
for the single baryon generated as an instanton. It is 
not yet known how to do this for dense baryonic matter 
although there has been an attempt toward that direc- 
tion [Ig . [l6ll . As the first step, we follow the strategy 
of Refs. [13, ^^ to integrate out all but the lowest vec- 
tor mesons p and w of the infinite tower of vector and 
axial vector mesons, making sure that hidden local sym- 
metry is preserved up to O(p^) in derivative expansion. 
The structure of single baryon has been worked out in 
Refs. [l|, Q using the hologra phi c QCD model of Sakai 
and Sugimoto (SS for short) [13, UM ^^^ ^-l^o with flat 
space in the SS model. This brings about a drastic simpli- 
fication from the infinite-tower theory but it has a great 
advantage: Once the pion decay constant /^ and the vec- 
tor meson mass mp are fixed from experiments, there are 
no undetermined parameters. As explained in Ref. [^, 
this comes about by the "master formula" which allows 
us to connect all the low-energy constants (LECs) of the 
(chiral) Lagrangian to the two constants /^r and nip. 

The most important feature that is found in Refs. ^^ 
was that the 0{p'^) terms in the Lagrangian, in particular 
the hWZ terms that carry information on the a; meson, 
could not be approximated by only a few terms - not to 
mention ignored - for the nucleon structure. This obser- 
vation suggested that the role of vector mesons, particu- 
larly, the w, would be crucial in dense baryonic matter. 
Indeed this will be what we will find in the present workjfl 



The objective of this paper is to put the skyrmions 
constructed as described in Refs. [l|, ^ on the crystal 
lattice as was done for the Skyrme model. Before stat- 
ing our principal result, we should mention a possible 
shortcoming in our model. When an HLS Lagrangian 
is implemented with baryons (that we shall refer to as 
BHLS), the theory in principle should give at the mean 
field approximation a relativistic mean-field model that 
describes nuclear matter properties as in the Walecka 
mean-field model. However as is well-known, it is es- 
sential that there be a scalar meson in the theory giving 
rise to the necessary attraction for binding [T^. (See 
Ref. [20| and references given therein.) Such a scalar, 
identified in Ref. [201 ^^ "dilaton," is missing in the HLS 
framework used here. How to introduce a scalar in the 
theory is highly problematic so we arc ignoring it in this 
paper but it may not be justified as we know from the 
original Skyrmc model that a correct nucleon mass can 
be obtained only when the Casimir contribution - which 
is 0{Nc) and captured in chiral perturbation theory at 
higher loop orders in the scalar channel - is taken into 
account. As we will discuss, some of the attraction is gen- 
erated by the w-p coupling in the crystal. However the 
missing dilaton degree of freedom may still be needed for 
more reliable calculations. We will return to this matter 
in the last section. 

The principal finding in this work is that the effective 
pion decay constant in medium /* decreases smoothly up 
to the density n^,^ at which a skyrmion fractionizes into 
two half-skyrmions and then stays non-zero, more or less 
constant up to the chiral restoration density n^. The no- 
table consequence is that the effective nucleon mass mj^ 
also stops scaling for n > n^,^- This is consistent with the 
large Nc consideration. In the large Nc limit, a nucleon 
in medium as a skyrmion scales as /* and hence remains 
more or less constant up to n^. What this means is that 
the nucleon mass could have a component that remains 
non-vanishing up to the chiral transition, which is rem- 
iniscent of the parity doublet picture of the baryon [21| 
with the nucleon mass given by 



iTT-N ^fna + A{{qq)), 



(1) 



where ttiq is the chiral invariant mass and A is the part of 
the mass that vanishes as (qq) — ^ (in the chiral limit) for 
n ^ Uc- The same result was obtained in a Rcnormal- 
ization Group (RG) analysis of hidden local symmetry 
Lagrangian with baryons in Ref. [20| and also in a phc- 
nomcnological study of the equation of state (EoS) for 
dense compact-star matter in Ref. [2a |. 

That rrig ^ implies that there is a part of the nucleon 
mass that is not generated by spontaneous breaking of 
chiral symmetry. The question then would be whether 
TOq reflects a fundamental feature of QCD as suggested in 



^ Note that the curvature is fixed by the gravity in the latter while 
undetermined in the former. As will be specified, we will deal 
with a "master formula" for any curvature. 

^ As will be shown explicitly, the minimal version of the HLS model 



HLSminCf, P, t^) exemplifies that it is a poor approximation to 
pick only one term out of three that are in the hWZ action. 



Refs. [23|, |2^ or an emergent symmetry via correlations 
in medium as in condensed matter as indicated in this 
crystal calculation. This raises the issue of the origin of 
the nucleon mass as we shall discuss below. 



with the chiral fields ^^ and ^^, which arc expressed in 
the unitary gauge as 



d = e 



B. 



o*'r/2/. 



e 



with 



TV ■ T, 



(3) 



This paper is organized as follows. In Sec. HI] we briefly 
sketch the derivation of the HLS Lagrangian from holo- 
graphic models that we are using in the exploration of 
dense matter. In Sec. lIIII we first discuss the main consid- 
erations in the investigation of the dense skyrmion matter 
using the face-centered cubic (FCC) crystal in the frame- 
work of HLS, then provide our numerical results of dense 
skyrmion matter and discuss their physical implications. 
We study the in-medium properties of hadrons with FCC 
crystal as the background in Sec. lIVI In Sec.|Vl we match 
our results to those obtained in (1) nuclear effective field 
theory (EFT) focused on the EoS for compact star matter 
and (2) RG-analysis of HLS theory with baryons. Fur- 
ther remarks relevant to the physics of rare isotope beam 
(RIB) facilities and conclusions are given in Sec. I VII 



where r's are the Pauli matrices. The covariant deriva- 
tive is defined as 



Df^^RX = (^M - «^p)Cfl,L 



(4) 



with Vf^ being the gauge boson of the HLS. This is the 
way to introduce vector mesons in the HLS, where the 



vector meson field V^ is 



with 



^P = ^ (^M + Pm) 



Pm = Pm 



\V2p- -pI 



(5) 



(6) 



Up to O(p^) including the hWZ terms, the most gen- 
eral HLS Lagrangian can be expressed as 



C 



HLS 



C 



(2) 



c 



(4) 



c. 



ij) 



with 



II. HIDDEN LOCAL SYMMETRY AND SINGLE 
SKYRMION SOLUTION 



C 



(2) 



fl Tr (a_L^a'[ 
1 
25 



«/^ Tr (aiipaji' 



;^Tr {V^,Vn , 



(8) 



Wc briefly describe the HLS Lagrangian that we shall 
use in exploring dense matter. The full symmetry 
group of the effective Lagrangian is Gfuii = [SU(2)^ x 
SU(2)^]ehirai X [U(2)]hls- Hcrc, [U(2)]hls is adopted as 
the "hidden local symmetry" to incorporate the p and 
u) vector mesonsQ The building blocks of the HLS La- 
grangian are the two 1-forms d||^ and dj_p defined by 



h^^ = Y^{D^S,ReR + D^i^eL). 



TSD,UeR-D,£.L£,l), 



2i 



(2) 



where /^ is the pion decay constant, a is the parameter 
of the HLS, g is the vector meson coupling constant, and 
the field-strength tensor of the vector meson is 



1^^. = a^K-a,v;,-i[i/^,K] 



(9) 



In the most general form of the O(p^) Lagrangian there 
are several terms that include two traces in the flavor 
space such as the jJiQ-yig terms listed in Ref. @ These 
terms are suppressed by Nc compared to the other terms 
in the Lagrangian and are not considered in the present 
work. Then the O(p^) Lagrangian which we study in this 
paper is given by 



£ 



(4) 



-(4)^ 



c 



(4)2 



(10) 



wher 



£ 



(4)y 



2/1 Tr 
+ y5Tr 
+ %{Tr 



t±„Lt± 



y2Tr 



^±m"-Li/"±"± 



"_L,,"'fd^„d|| 



%Tr 

+ Tr 



"_LA'"j-^"||"li 



"_LM"lk"j^"|| 



ysTr 
l+2/9Tr 






y4Tr 



"llM"lk"ll"li 



''^ti'^^y'- 



"iM^ik^'i^ii 



(11) 



4 Here for simplicity we are assuming that flavor (7(2) symmetry 
applies to the p and ui in medium as it seems to do in free space. 
It is found by the RG analysis in Ref. [20l | that the U{2) must be 
broken down appreciably in medium, so there may be a caveat 



in what wc obtain here. 
^ Another example of this kind is Tr[Q{^]Tr[Q:||p] that generates 
the mass difference between the p and oj mesons. 



C-(i)z = «^4Tr 



V^ya'\_a''^ 



iZgTr 



Vu^dlfaj; 



(12) 



Finally, the anomalous parity liWZ terms are written as 



£« 



167r2 ^ 



^i -^ i 5 



where 



£i = i Tr [aLttR - a%^a^ , 
£2 = i Tr [ttLaR^LaR] 7 
£3 = Tr [Fy {ai^a^ ~ a^a^^)] , 
in the 1-forni notations with 



(13) 



(14a) 
(14b) 
(14c) 



a, 



a„ 



a 



'■L — "II ~ "±i 

i^y=dT/-iT/^ (15) 

As for the low energy constants F^, a, (7, j/i, Zi, and Ci, 
we use the general "master formula" proposed in Refs. [l|, 
[2| from a class of holographic QCD models, namely, 

5 

2/1 = -2/2 = -^cGym ((1 + V^i - V^o)') 

2/3 = -2/4 = -^cGym (v^? (1 + ^if 

2/5 = 2yg = -2/g = -2A^cGym (i^Ii'o) : 

2/6 = - (2/5 + 2/7) > 

2/7 = 27VcGym (V-i (1 + V-i) (1 + ^1 - V-o)) , 

Z4 = 27V,Gym(V'i(1 + V'i-^o)), 

Z5 = -2iV,GYM(V'?(l + V^l)), 

1 
2 



V-o^i 0^0 + «^i 






6 



1,21,2!, 1 
-2^0 + 6^? +2^^1 + 2 



(16) 



where the smallest nonzero eigenvalue X^ and its corre- 
sponding eigenfunction ipi{z) satisfy the eigenvalue equa- 
tion, 

- Kr\z)d. [K2iz)d.4'niz)] = A„^„(z) (17) 

with Ki{z) and K2{z) being the warping factors of the 
5-dimensional space-time. In Eq. (J16p we have made the 
following definitions: 

/•OO 

{A) = / dzKi{z)A{z), 



m) 



dzA{z) 



(18) 



In the present work, we take as in Refs. [ll, l2| Ki{z) = 
K-^/^{z) and K2{z) = K{z) with K{z) = 1 + z^ corre- 
sponding to the Sakai-Sugimoto model [10, llll • There- 
fore, once the hQCD parameters Mkk and Gym are de- 
termined by using the real world parameters, for exam- 
ple, such as /^ and ijip, all the rest parameters of the 
effective Lagrangian except the parameter a are fixed 
through the master formula. In this work, we will work 
with the parameter set determined by the epirical values 



U = 92.4 MeV, 



rUn 



775.5 MeV. 



(19) 



Note that, due to the relations between y^'s, the to 
meson can couple to the others only through the hWZ 
terms. We will denote the w-decoupled HLS Lagrangian 
with Cj = for all i as HLS(7r, p) in order to distinguish 
it from the full Lagrangian IILS(7r, p, uj). If we keep only 
the terms up to 0{p^) and only one term among the hWZ 
action that results if one takes the equation of motion for 
the p meson and the large mp limit in the hWZ action in 
Eq. ([7]o and use the same values for the parameters, the 
Lagrangian is reduced to the "minimum model" adopted 
in Ref. [25|, which is denoted as IILSmin(7r, p, w). 

In Refs. [l|, H, the single skyrmion solution was found 
under the hedgehog configuration for various parameter 
sets. For HLS(7r, p, oj) with the parameters determined 
by Eq. ([H]), the soliton mass Msoi is 1184 MeV and 



the root mean square winding number radius \/ {r'^)w 



is 0.433 fm 
yields M^oX 



Compared to the HLSn 



{.r'^)w 



i{tt,p,uj) which 
= 0.540 fm, we 



1407 MeV and V> 
can see that the 0{p*) terms provide about 300 MeV re- 
duction in the soliton mass and about 0.1 fm reduction 
in the size. This means that the n-p interactions incorpo- 
rated through the O(p^) terms and the ir-p-uj interactions 
through the hWZ terms provide such an effective attrac- 
tion. As we will see later, the reduction in the size of the 
skyrmion plays the most important role in determining 
the critical density n^,^, where the half-skyrmion phase 
transition happens. 

On the other hand, if the w is decoupled, Mgoi is re- 
duce d to 83 4 MeV and the skyrmion has a very small size 
as \/{r'^)w — 0.247 fm. This emphasizes the role of the 
repulsive u meson in the HLS Lagrangian. 



As mentioned, this limit is at odds with the basic premise of 
HLS il. 



TABLE I. Symmetries of the FCC skyrmion crystal. 



Reflection 
(yz-plane) 



3- fold axis 
rotation 



4-fold axis 
(z-axis) rotation 



Translation 



{x,y,z) 

U — a + iT ■ TT 

Pi = SaipPp 



{-x,y,z) {y,z,x) 

(a, — 7ri,7r2,7r3) (cr, 7r2, tts, tti) 

{-P1,P2,P3) {p2,p3,pl) 

UJO UJO 



{x,z,-y) 
{a, TTijTTa, — 7r2) 

(P1,P3,-P2) 
UJO 



{x + L,y + L,z) 
{a, — TTi, — 7r2,7r3) 

(-P1,-P2,P3) 
UJO 



III. SKYRMION MATTER FROM THE 
HIDDEN LOCAL SYMMETRY 



Since the skyrmions are classical solitons, they pro- 
vide a natural framework to study dense matter. We 
can make the dense skyrmion matter by putting more 
and more skyrmions into the space without worrying 
about such a serious problem as the sign problem in 
the lattice QCD. What we consider is a crystal made 
of skyrmions. In the lowest energy configuration, the 
skyrmion at each lattice site should be arranged in such 
a way that the nearest skyrmions have the maximum at- 
traction. In the study with the original Skyrme model, it 
turns out that two skyrmions apart from each other has 
the maximum attraction when one of them is rotated in 
the isospin space with respect to the other by an angle 
TT about the axis perpendicular to the line joining them. 
In Rcf . [2 a ] , it was shown in a simple approximation that 
such a configuration is also most attractive in the model 
of HLSniin(7'', p, Lo). If wc assumc that it holds in the case 
of HLS(7r, p) and HLS(7r, p, w), then the FCC crystal con- 
figuration for the TT, p, and u should have the symmetries 
with respect to the reflection and translations as summa- 
rized in Table HI 

The classical solutions for n, p, and uj mesons satisfy- 
ing the symmetries for the maximum attraction and the 
constraints for each FCC box to carry a given baryon 
number can be obtained by ri sing the Fourier expansion 
method developed in Refs. [2a,[23 for the original Skyrme 
model with the pious only, then generalized in Ref. [23 
for the model with vector mesons. The main idea is to 
introduce the un- normalized fields 4>f/ {M = tt, p), each 
of which can be expanded in terms of the Fourier series 
as 



^o"^ 



ir 



/^ Pabc cos{7rax/L) cos{nby / L) cos{'kcz/ L), 

ahc 

y ai^i^i sm{TThx/L) cos{TTky/L) cos{nlz/L), 

hkl 



"2 ^/ cthkiCOs{TTlx/L)sm{TThy/L)cos{TTkz/L), 

hkl 

63^ ~ y^ ctfiii cos{7rkx/L) cos{7rly / L) sm{7rhz/L), 

hkl 



and thereafter normalized as 



aJ>/ 



IM 



E,(^^)" 



(20) 
(21) 



for each M = tt, p. Here, 2L is the size of a single 
FCC box containing four skyrmions inside. Thus, the 
skyrmion number density n is given by n = 1/{2L^). If 
we take skyrmions as nucleons, the normal nuclear mat- 
ter density Uq = 0.17/fm'^ corresponds to L ^ 1.43 fm. 

Finally, U{= ^^) and pf can be expressed in terms of 
these normalized fields as 



U 



0^ 



9P7 



£abc<l>bdi<l>^ 



1 



(22) 



The isoscalar field ojq has the same symmetry proper- 
ties as the (/jq and therefore could be expanded as 

Wq = yj /3^^c cos{T:ax/L) cos{nhy/L) cos(7rcz/L)(23) 

ahc 

with the expansion coefficients /?J^f,c- 

The Fourier coefficients and their indices are subject 
to the topological structure and the constraints from the 
symmetry properties listed in Table HI From the sym- 
metry properties one can easily check that in the above 
equations (a, 6, c) should be all even or all odd integers, 
and (fc, /) should be all even (odd) if h is odd (even). 
Furthermore, to provide the correct topological structure 
(winding number) to the configuration, the expansion co- 
efficients /3J^^ must satisfy the constraint. 



even 



Pab. 



0. 



(24) 



As for the pion and p field configurations, all the Fourier coefficients are determined by minimizing the skyrmion 
energy per baryon, E/B, which, from the Lagrangian ([7]), is expressed as 



E/B = (E/B) 



0{p^ 



(E/B) 



0(pi 



(E/B) 



hWZ, 



(25) 



where 



d?i 



-f {a±, ■ d_L^) - -^uj^ + ^ (a||i • dni) ~ -d,ujod,ujo + -^V,j ■ V^ 



(26) 



{EjB) 



0(pi 



d^x 



(.^/-DJhWZ — T 



4 V 327r2 



4 



2/5 



{( 



q; I i X CKiu" • [OL \ j X an 



+ -J {oi^i X q:_lj) • [a.\\i X aiij) 



(^aj_i X apj • (^a^j x an^- 

z 



)} 



^Vy • (aj_, X aj_j) - ^Vy • (ay, x ay^-j 



d^x 



-Ci^qEijA; [a^i ■ (ajij x q;||j,) + a^i ■ {ot^^ x a^k)] 



,(27) 



(28) 



with the integration being taken over the space inside an 
FCC box of size L, and aj_i, q;||j, and Vij are expressed 



in terms oi (p'^'P as 



'±i 



^^r 



rd^^O 




l + '/'o 

(a.(/." X (j)'') {d,r X cpP) 



(1 + </>?) 



(1 + 0S) 



(29) 
(30) 



1 









Xdrr^d.rT 



1 



+ [(.9pOx(5P,)]° .(31) 



On the other hand, the simple minimization procedure 
of E/ B to find the expansion coefficients /3^^^ would lead 
to the trivial solution Wq = 0. This is because there are 
no constraints for the |iabc^^ ^^-d the a; always provides 
repulsive interactions. Therefore, the w expansion coef- 
ficients should be determined by satisfying the equation 
of motion for uj{= Wq/Ztt); viz. 

(-a,9, + ml) uj{r) = r- d,jr, (32) 

where J" and diJf coming from the hWZ terms are 






327r 



2^ijft 



(ci -t-C2)Q:±i • (q:||j X q;||;,) 
+ (ci - C2)q:±i • {ot^j X a^fc) 
- 2c3V;j • a_Lfc , 



5^c - - 

^ Co £,-,fc diu- • q: I I. . 



(33) 



I 

panded in terms of the Fourier series as 

J" — diJf = 2_\ "fabc cos(nax/L) cos{Trby/L) co&{-ncz/L). 

abc 

(34) 
Then, thanks to the linear nature of the equation of mo- 
tion for the a;, /S^^^ can be easily obtained as 



"fabc 



abc 



{Tr/L)^{a^ + b^ + c^)+ml 



(35) 



Thus, in the numerical minimizing processes only the 



''hkl 



and 



abc 



are taken as adjustable parameters. 



Shown in Fig. [T] are the numerical results for the de- 
pendence of E/B and (a) on the FCC size parameter 
L. Here, (a) is the average of the chiral field U over the 
space that a single skyrmion occupies and it vanishes in 
the half-skyrmion phase. Here, we present the numerical 
results for HLS(7r, p, w), HLS(7r,/ci) and HLSmin(7r, p, w) 
by solid, dashed, and dash-dotted lines, respectively. The 
position of the normal nuclear matter density is denoted 
by a vertical dotted line. 

One can see that in HLS(7r,p, w), the half-skyrmion 
phase transition takes place close to the normal nuclear 
matter density Uq. Comparing to that of the minimal 
model HLSmin(7'',/c, w), we see that the inclusion of the 
O(p^) interaction terms and the other w-p-uj interactions 
through the hWZ terms makes n^,2 larger by a factor 
of 1.7, which is quite a noticeable effect. This higher 
value for n^,^ in HLS(7r,p, w) compared to that of the 
minimal model may come from the fact that the size of 
the single skyrmion is smaller in the former and that the 
additional interactions in HLS(7r, p, uj) weaken the repul- 
sive interactions from the uj. This analysis is supported 
by the results from HLS(7r, p), where the half-skyrmion 



phase transition occurs at L ^ 0.8 fm, i.e., n^ 



/2 



6 no. 



Due to the crystal structure, the source term can be ex- much higher than with the uj. This also accounts that 
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FIG. 1. E/B and (a) of the minimum energy configuration 
for a given L. Here and in all the figures that follow, L de- 
creases to the right indicating increase in density. The numer- 
ical results for HLS(7r,p,tj), HLS(7r,p), and HLSmin(7r,p,a;) 
are presented by solid, dashed, and dash-dotted lines, respec- 
tively. 



in HLS(7r, p), the absence of the lo reduces the skyrmion 
size to nearly half of that of HLS(7r, p, lo). 

The density n^^j^ at which E/B is the minimum turns 
out to be larger than n^,^ for all three models con- 
sidered in the present work. As one can see from 



/2 



Fig. [T] the minimum is achieved at n^^:^^/n-^ 
HLS(7r,p, oj) and in HLSniin(7i', P, w), but at ?^min/"'i/2 '" 
1.5 in HLS(7r,p). The binding energy per baryon at this 
minimum energy configuration can be also read from 
Fig. m as - 150 MeV, - 100 MeV, and - 50 MeV 
for HLSinin(7'',p, w), HLS(7r, /9, w) , and HLS(7r,p), respec- 
tively. As a complete theory of nuclear matter, the mini- 
mum at rijfjin should represent the nuclear matter ground 
state. That n^:^^ > rig and the binding energy at n^^^ is 
larger than the empirical value ^ 0.17 fm~'^ signals that 
there is something missing in the present skyrmion crys- 
tal description for the EoS of nuclear matter. We will 
discuss the plausible causes for this defect and possible 
remedies in the last section. 

It is of theoretical, though perhaps academic, inter- 
est that the E/B versus density of HLS(7r,/9) is (nearly) 
flat. One can easily understand this feature, which arises 
from the important role of the vector mesons in nu- 
cleon structure. In Refs. [2^, [2^, the 5-dimcnsional BPS 
theory - that can be considered as the Sakai-Sugimoto 
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FIG. 2. Breakdown of E/B. The solid, dashed, and dash- 
dotted lines represent {E/B)Q(^p2)^p, {E/B)Q^pi^^p, and 
{E/B)ui, respectively. 



hQCD model with flat metric - is transformed to a 4- 
dimensional infinite-tower meson theory and the instan- 
ton in that theory that satisfies the Bogomolnyi bound 
condition is then transformed to the skyrmion. It is 
shown that, when only the pions are considered, the 
skyrmion carries an excess energy of about 25% relative 
to the Bogomolnyi bound. However, when the lowest 
SU{2) vector meson, the p meson, is included, the en- 
ergy is drastically reduced, arriving almost at the Bo- 
gomolnyi bound. The next vector meson, which is the 
axial ai, brings it even closer to the bound. What is 
remarkable - and highly relevant to our work here - is 
that the lowest vector meson p alone brings the energy 
nearly to the bound. This is highly reminiscent of the 
vector dominance of the nucleon isovector electromag- 
netic form factor by the p (with a small correction from 
the next vector meson p') p^ - [l4l . [3g |. Since the HLS 
Lagrangian used in our work is constructed from the 5D 
Sakai-Sugimoto model in a similar way, the difference be- 
ing the background curvature, the inclusion of the p, and 
the higher order terms in the Lagrangian makes the single 
skyrmion mass almost saturate the Bogomolnyi bound. 
Consequently, there is very little interaction energy left 
in the skyrmion made of the pion and the p to disturb 
the energy versus density curve. 

We should stress that the approach to the Bogomolnyi 
bound by the p - and higher isovector vector mesons is 
obstructed by the t/(l) vector meson w (and its tower). 
It is this uj that plays a crucial role in nuclear dynamics 
as will be explained below. 

As emphasized and proved in Refs. [l|, 0], our calcu- 
lation of skyrmion properties is independent of the HLS 
parameter a\j The same argument can be made here, 



'^ In standard HLS, a is a free parameter which is normally taken 
to be 1 < 1 !5 2 [aigl. In free space a ~ 2 is preferred, but 
in a hadronic medium at high temperature and/or density, one 



i.e., E/B and the expectation value (ct) are independent 
of the choice on the parameter a. Again, as in Refs. [l|,yl, 
we could confirm the a-independence of E/B and (a) in 
the present work. 

Figure[2]presents the contribution of each term in E/B 
given in Eq. ([25]). Here, {E/B)o{p2)^^p, {E/B)o{p4)^,p, 
and [E/B)^ are given by the solid, dashed, and dash- 
dotted lines, respectively. First of all, one can find that 
{E/B)o{pi)TT,p is negative for all density range involved, 
showing clearly that the account of the entire 0{p^) terms 
and vector mesons brings in a sizable attraction, a feature 
that has not been exposed before in the literature. What 
is more striking is its weak density dependence compared 
to the other two contributions. Next, {E/B)o{p2)Tr,p and 
(E/ B)uj show distinctly different density dependence in 
the skyrmion phase and in the half-skyrmion phase. In 
the single skyrmion phase, {E/ B)o(p2)TT.p dominates over 
{E/B)^, with their ratio [E / B)o[p2)T:^p/ {E / B)^ even in- 
creasing from ^ 2 to '--^ 3 as the density increases whereas 
in the half-skyrmion phase, [E / B)o{p'2)tt,p starts to de- 
crease while {E/B)^ starts to increase. Then, at some 
density the repulsive interaction due to the w becomes 
dominant. 



IV. HADRON PROPERTIES IN THE FCC 
CRYSTAL BACKGROUND 

A. In-medium modification of meson properties 

As stressed, the power of the skyrmion-mattcr model 
is that it describes baryons and mesons in a unified way. 
This, therefore, makes it a natural framework to in- 
vestigate the in-medium modification of meson proper- 
ties [J, [31| . The Lagrangian ([7]) itself could describe the 
meson dynamics in baryon free space. If we expand the 
Lagrangian in terms of the fluctuating meson fields about 



their vacuum values, that is, [/ = 1, p° 
Lagrangian ([7]) becomes 



0, the 



-m»-'' 



interaction terms, 



(36) 



m?p = w},,. 



where we have used ag^f^ 

On the other hand, what we have obtained in the pre- 
vious section is the classical solutions for U, p°, and cj^, 
which are the minimum energy configurations for a given 



baryon number density. If we incorporate the fluctua- 
tions with respect to these classical solutions, then the 
fluctuating fields describe the corresponding mesons and 
their dynamics in dense baryonic matter, which is the 
main idea of Refs. [3|, [3l[. We denote the minimum en- 
ergy solutions as C/(o) = C(^o)L^(o)it;' /^?' > ^'^^ ^1°^ ^^'^ 
introduce the fluctuating fields on top of the classical so- 
lutions as 



^L,R — C(o)L.i?^L.fl 



Wu 



Hp 



>' 



(37) 



where ^\ = ^r ^ l = cxp(zTa7ra//^), p^, and ibp stand 
for the corresponding fluctuating fields. By substitut- 
ing these fields into the HLS Lagrangian and taking the 
space average for the background fields configurations as 
denoted by (•), we are led to 



C 



{a- 



1 
1 

V 
1 

4 



„ 1 



'^h) 



9^7r"5^7r" 



i.^ypP''^ 



p^^f"" + -mtujpuj^' + 



(38) 



where we have considered the 0(p^) terms only. The 
extra factor in front of the pion kinetic term can be ab- 
sorbed into the redefinition of the pion fields as 



C = exp(iTQ^*//;), 



(39) 



with 






l-(«-l)3(l-{-fo)))- (40) 



Since we have considered only the 0{p^) terms, there 
does not appear any modification factor for the p meson 
and f^/f-K becomes dependent on the parameter a. If 
we include O(p^) terms and one-loop corrections from 
the 0{p^) terms correctly, we believe that such an a- 
dependence will disappear or be weakened. 

In Fig. [3l we present /^//tt as a function of the FCC 
box size, where we have fixed the parameter a to 2|f| 
Again, the density dependence of f^/ f-r: shows the dif- 
ferent behavior in the single skyrmion phase and in the 
half-skyrmion phase. In the single skyrmion phase, /^//tt 
decreases as the baryon number density increases. On the 
other hand, in the half-skyrmion phase, it stays in a non- 
vanishing value around 0.65. Due to this non-vanishing 
f* , the half-skyrmion phase was interpreted as a kind of 
pseudo-gap phase in Ref. [S^ ■ 



gets a ~ 1 [3|. That the parameter a does not figure in physical 
observables in this theory raises the question as to how the vector 
dominance, a powerful notion in hadron physics, manifests itself 
in the formalism, in particular in medium. 



We include only the 0{p^) terms, which can reasonably repro- 
duce the meson properties only when we take a = 2. 
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FIG. 3. Naive estimates of /jj//,r as a function of crystal 
size with the skymiion matter background. Notations are the 
same as in Fig. [l] 



B. In- medium modification of baryon properties 

In our HLS Lagrangian all the parameters are deter- 
mined by the master formula (J16p given in the matter-free 
space. As the simplest approximation in applying it to 
dense matter, we shall assume that the in-medium modi- 
fication of all the parameters involved, such as the 0(p'') 
terms, can be calculated in terms of the in-medium con- 
stants f* and m*. We denote them with asterisk. Now 
given these starred parameters, we can solve for the single 
skyrmion solution as in Ref. [^ to obtain the in-medium 
modification of the baryon mas^j in dense medium. To 
get the leading corrections, it suffices to use the O(p^) 
Lagrangian with a = 2 as suggested by phenomenology 
to calculate the two in-medium constants, in which case 
the most important effect is lodged in the pion-decay con- 
stant. The effects of the 0{p*) terms so calculated can 
then be straightforwardly incorporated with the starred 
quantities. In principle, we should rewrite the O(p^) 
terms in medium to obtain a relation which may be dif- 
ferent from the free-space master formula, but this will be 
complicated considering that Lorentz invariance is bro- 
ken in medium. Intuitively, however, we do not expect 
the result to be significantly modified. 

Shown in Fig. |4] is the in-medium single skyrmion 
masses calculated in the approximation described above 
as a function of the FFC size parameter L. We also 
present the A-A^ mass difference in Fig. [5] and the soliton 
root-mean-square (rms) radii in Fig. [6] 

Note that the behavior of the crystal size dependence 
of the soliton mass shown in Fig. |4] is the same as the 
behavior of f* shown in Fig. [3l This can be understood 
by that in the skyrmion model the soliton mass scales 



^ In the previous section, we have studied the skyrmion energy per 
baryon E/B. This does not represent the quasiparticle (baryon) 
mass since it includes the interactions among the baryons. 
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FIG. 4. Crystal size dependence of soliton mass with the FCC 
crystal background. Notations are the same as in Fig. [T] 
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the FCC crystal background. Notations are the same as in 

Fig.m 



roughly as Afgoi '^^ /^, which follows simply from the 
scaling of the nucleon mass in the large Nc limit as in 
matter-free space. Figure [5] shows that the crystal size 
dependence of the A-A^ difference Aj\/ changes a lot de- 
pending on the existence of the to meson. In the calcula- 
tion with the uj meson, Am decreases up to n-^ /j beyond 
which Am becomes nearly constant. However, in the 
case of HLS(7r, p), Am is unsealing up to n-^^,^ but, af- 
ter n^,2, A A/ increases dramatically. Our finding that 
the nucleon mass, which decreases with increasing den- 
sity in the skyrmion phase, stops dropping at n^,^ and 
stays constant in the half-skyrmion phase is the most 
significant result of this paper and, up to our knowledge, 
has not been obtained in other models. In the case of 
HLS(7r, p, oj), the in-medium soliton mass approaches to 
~ O.GMsoi. As will be discussed in the next section, 
this finding will get some support from other approaches, 
both theoretical and phenomenological, and implies con- 
sequences on upcoming experiments. 

Similarly, the skyrmion radius is strongly infiuenced by 
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FIG. 6. Crystal size dependence of soliton radii with the 
FCC crystal background. The upper panel is for the winding 
number rms radius and the lower panel is for the energy rms 
radius. Notations are the same as in Fig. [T] 



the presence of the ut meson as seen in Fig. [HI In the case 
of HLS(7r,p), the skyrmion radius is small and constant 
at varying density. This is simply because in the absence 
of the a;, the skyrmion is almost point-like and remains 
stable with the crystal size. This feature was already 
obvious with the instantons in the hQCD model. 

It is noteworthy that, while the impacts of the 
skyrmion-half-skyrmion transition can be seen in all 
three models for the masses and splittings indicating the 
effect of topology change in medium, it does not affect 
the skyrmion size when the u is absent, which is in a 
stark contrast to the case when the lu is present. 



V. ORIGIN OF NUCLEON MASS 

The principal finding in this paper, as mentioned be- 
fore, is that simulated on the crystal lattice, the nucleon 
mass decreases smoothly as density increases up to n-^,2 
at which skyrmions fr actionize into half-skyrmions and, 
in the half-skyrmion phase, the dropping of the nucleon 
mass stops and the mass remains constant going toward 
to the chiral restoration. In the skyrmion phase with 



n < n 



1/2' 



chiral symmetry is broken with the order 
parameter (qq) ^ 0. Hadrons are massive apart from 
the Nambu-Goldstone bosons, the pious. In the half- 
skyrmion phase, (qq) = on the average, but chiral sym- 



metry is not restored since hadrons are still massive and 
there exist pions. There seems to be no obvious order pa- 
rameter characterizing this state apart from the presence 
of the half-skyrmion structure, the changeover involving 
a sort of topology change. 

Apart from the changes in the baryon properties dis- 
cussed above, there are several other remarkable happen- 
ings at n^,^. For instance, in Ref. j33| was found that an 
antikaon K~ propagating in the skyrmion matter under- 
goes propitious drop in its effective mass at n^,.^, which 
may have an influence on kaon condensation in compact- 
star matter. Another phenomenon concerns the symme- 



try energy of nuclear matter. At 



''1/2- 



the symmetry 



energy is found to develop a cusp structure which leads 
to drastic modification of the nuclear tensor forces [3J| 
that contribute importantly to the structure of asymmet- 
ric nuclei and large mass neutron stars ^^. We believe 
that all these phenomena are closely linked together and 
indicate the possible source of the proton mass. 

There are several indications that the nucleon mass 
does stay constant going toward the chiral transition 
point. It points to a possible origin of the proton mass. 
The first indication comes from a phenomenological ap- 
plication of the skyrmion-half-skyrmion transition to the 
EoS of nuclear and neutron matter j22| . The idea is to 
translate what is observed in the skyrmion description 
to the density dependence of the parameters of the HLS 
Lagrangian, namely the masses and coupling constants 
scaling with density, with which a nuclear effective field 
theory is constructed. The topology change at n-^^,^ im- 
pacts on the changes in the parameters of the EFT La- 
grangian and the drastic effects mentioned above reflect 
on nuclear forces, in particular the tensor forces. It is ob- 
served in the EoS at high density that the nucleon mass 
should stabilize at n^,^ to m*pf ~ O.Smj^ while the ef- 
fective ujNN coupling g* jvat can drop but very slowly in 
density, constrained by the constant nucleon mass. Na- 
ture seems to be exposing an intricate interplay between 
TO^ and the coupling .g*^jv or more generally the w de- 
gree of freedom. 

Another indication comes from an RG analysis of nu- 
cleon properties in HLS with baryons, i.e., BHLS, as scale 
is changed. It is found in Ref. [20| that as one approaches 
the dilaton limit fixed point at which the nonlinear sigma- 
model phrased in HLS transforms to a linear sigma model 
of Gcll-Mann-Levy type, the nucleon mass that falls in 
dilute environment stops falling at n^ and then stays 
constant afterwards, strongly correlated with the wNN 
coupling. It is suggestive that n^ corresponds to n^,^ m 
the crystal description. 

That the nucleon mass remains non- vanishing as one 
approaches the chiral transition point can be captured 
by the nucleon parity-doublet model [2l[ with the two- 
component nucleon mass ([1]) with a substantial nin ~ 
(0.6 — 0.8) 771 jv that remains non- vanishing. In Ref. j20| . 
this model was used to show that the RG fiow to the 
infrared fixed point, i.e. dilaton-limit fixed point, does 
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indeed indieate that there is a changeover at some den- 
sity for the nucleon mass in a close interplay with the 
w-nuclcon coupling. The non-zero TOq can be interpreted 
as reflecting a symmetry intrinsic in QCD as suggested 
by Glozman et al. [2j, [2J| . This interpretation suggests 
a mechanism for a bulk of the mass generated differently 
from the spontaneous breaking of chiral symmetry, pre- 
sumably tied to chiral invariant gluon interactions. The 
interesting question is whether one cannot identify the 
agent for the non-zero uiq with one that triggers the 
spontaneous breaking of scale symmetry that makes the 
dilaton condense. 

On the contrary, one can consider mg representing a 
symmetry that is "emerging" from strong correlations 
in nuclear interactions, which has nothing to do with a 
fundamental symmetry of QCD. This is what is coming 
out of our skyrmion crystal picture. In either case, it is 
an interesting issue in nuclear physics, intimately linked 
to the origin of the nucleon mass. 



VI. REMARKS AND DISCUSSIONS 

In this paper, to investigate dense baryonic matter, we 
simulated on crystal lattice the HLS Lagrangian written 
to 0{p'^) in the derivative expansion with all the param- 
eters of the Lagrangian fixed except for the pion decay 
constant f^^ and the p meson mass nip. This comes about 
in integrating out all vector mesons in the infinite tower of 
vector and axial vector mesons from a generic 5-D hQCD 
action in a curved space plus the Chern-Simons term, 
leaving for relevant degrees of freedom only the lowest 
vector mesons p and w in a hidden local symmetric form, 
obtaining a master formula that relates all 0{p'^) low en- 
ergy constants to two constants that figure in hQCD La- 
grangian, i.e., /^ and mp. The dependence on the famous 
a parameter which figures importantly in the usual HLS 
theory gauge equivalent to non-linear sigma model, such 
as the KSRF-type relation niy — ag^f^, can be elimi- 
nated in a way analogous to a gauge parameter. Physical 
quantities are independent of a. 

The most notable findings in this work are (1) the 
prominent role of the cj meson encapsulated in the homo- 
geneous Wess-Zumino term (equivalently the U{1) gauge 
field in the Chern-Simons topological term in hQCD) in 
the structure of both elementary nucleon and dense bary- 
onic matter and (2) the presence of topological change 
from skyrmions to half-skyrmions at a density relevant 
slightly above normal nuclear matter density and higher 
density. 

The first which plays a crucial role in giving correct 
structure to the nucleon controls the location of the 
topology change n-^^,^ and accounts for the large constant 
mass to the nucleon that remains non-zero at the chiral 
restoration which is locked to the pion decay constant 
effective in the half-skyrmion phase. This feature com- 
bined with in-medium properties of mesons introduces a 
dramatic change in the nuclear tensor forces at a density 



above nuclear matter and gives rise to a stiff EoS that 
could accommodate large neutron star masses such as the 
1.97Mq object found recently [22|. Highly relevant to the 
forthcoming RIB machines is that the drastic change in 
the nuclear tensor forces caused by the topology change 
we found at a density near that of nuclear matter can 
have a strong consequence on the shell evolution of asym- 
metric nuclei. It has been shown that the monopole ma- 
trix element of the tensor forces follows closely the shell 
evolution of single-particle states [331 ■ Thus any changes 
in the tensor forces could be probed by looking at the 
shell evolution. 

We should point out a few caveats and approximations 
made in the paper that require further work. The most 
glaring defect of this work is that the minimum of the 
energy density of the system occurs at too high a density 
and with too large a binding energy. One reason may 
be that our approximations are anchored on large Nc 
considerations and 1/Nc corrections cannot be ignored 
in medium. This is a general problem with skyrmion 
physics which relies on the Lagrangian valid for large Nc- 

Perhaps more important is the mechanism that leads 
to a correct equilibrium density and binding in EFT ap- 
proaches to nuclear matter such as RG ViowK [Sa]- There 
one finds that without irreducible repulsive three-body 
forces, the system saturates at much higher density and 
with much larger binding energy, very similar to what 
is found in our theory. It is possible that the skyrmion 
crystal approach, while containing reducible many-body 
effects, misses a mechanism that corresponds to this ir- 
reducible 3-body force effect. 

Now one can ask whether the predictions made in a 
theory that fails to reproduce the correct properties of 
normal nuclear matter can be trusted when applied to 
higher density. There is no clear answer to this question 
for the moment. However it may be that the quantity 
that depends on topology is robust. It is seen in the 
calculation that although the location of ri-^^,^ depends 
on dynamical details, the existence of the skyrmion-half- 
skyrmion transition itself is generic. The RG analysis of 
Ref. (201 provides an independent support to this result. 
The issue then is where the transition takes place, and 
this can be ultimately checked by experiments. There is 
an indication that n^,^ is near Uq and the phenomenology 



presented in Ref. 



hints that n 



1/2 



(1.5 -2.0) no is 



compatible with nature, a region that can be accessed by 
the future RIB machines and is also highly relevant to 
compact stars. 

As already mentioned, one important missing ingredi- 
ent is scalar degree(s) of freedom in the Lagrangian. In 
the Skyrme model (with pious only), the Casimir con- 
tribution to the nucleon mass which comes at 0{N^) 
is attractive and significant, of order of ~ —0.5 GeV. 
This cannot be calculated accurately, so is generally ig- 
nored in the community. In fact, it is the same 0{N^) 
order effect that gives rise to the binding of K'^ to 
the skyrmion that describes the hyperons in the Callan- 
Klebanov model [37|. It is very likely that this (or part 
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of this) effect can be captured by a low-mass scalar as 
in Walecka's mean field model for nuclear matter. We 
have seen in our model HLS(7r,/9, w) that the lo-p-i: cou- 
pling generates - via hWZ terms - a significant attrac- 
tion in the nucleon, but there could be further attrac- 
tion that gets enhanced in dense matter. How to im- 
plement scalars into HLS(7r, p, w) remains a problem to 
be worked out. A scalar in the form of dilaton associ- 
ated with spontaneously broken conformal symmetry as 
exploited in Ref. [20| could provide an approach to the 
problem. 

Among the approximations we made, we note that the 
crystal-size (or density) dependence of /^ is obtained 
with the O(p^) terms of HLS. A complete calculation 
of this quantity up to O(p^) should include not only the 
tree- level O(p^) contributions but also one- loop correc- 
tions with the 0{p^) terms. But the latter are suppressed 
by the Nc counting. Moreover, since we only consider 
fiuctuations from the 0{p^), both the p and w meson 
masses are the same as their vacuum values and indepen- 
dent of the density. In a complete calculation including 
the full O(p^) contributions which contains not only the 
tree- level O(p^) terms but also one-loop corrections from 
the 0{p^) terms, these masses will change with density, 
although the changes will be suppressed by 1/Nc- 

Similarly for the density dependence of the soliton 
properties, only the in- medium modified /* is considered 
since the other input parameter nip is the same as its 
vacuum value in the present approximation. Therefore, 
the scaling behaviors of the soliton properties agree with 
the large Nc argument. In a complete calculation up to 



O(p^) including the one-loop corrections from the O(p^) 
terms, these results might be changed but again they will 
be suppressed by the chiral counting or the N^ counting. 
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